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An experimental foundation for electromagnetic shower formation in the 
geomagnetic field 

U.I. Uggerh0j a 

a Department of Physics and Astronomy, University of Aarhus, Denmark 

At very high energies the Earth magnetic field, although in many other connections considered rather weak, 
acts as a strong field on photons and leptons. This paper discusses the intimate connection between this effect 
and the corresponding 'strong field effects' observed at accessible accelerator energies in aligned single crystals. 
As such, these effects constitute an experimental verification of the theoretical basis used for simulations of the 
development of electromagnetic showers in magnetic fields, in particular the geomagnetic field. A short discussion 
of more general aspects of the shower development in the fields present at different distance scales is included. 



1. Introduction 

With the advent of facilities like the Pierre 
Auger Observatory [1 for the detection of ultra 
high energy cosmic rays of energies in the EcV 
(10 18 eV) region and orders of magnitude above, 
pair production and photon emission in the mag- 
netic field of the Earth become increasingly rel- 
evant. These phenomena have been studied as 
early as the late thirties by Pomeranchuk > later 
by McBreen and Lambert in the early eighties 3 
and recently e.g. by Stanev and Vankov 0] as 
well as in extended air shower (EAS) simulations 
by Plyasheshnikov and Aharonian 5 . An im- 
portant issue in this context is the possibility of 
distinguishing photon-initiated EASs from those 
initiated by protons or heavy nuclei as this dis- 
tinction may shed light on the question of 'top- 
down' (topological defects, massive X-particles) 
or 'bottom-up' (acceleration of known particles) 
mechanisms |6I7| . Another issue is the directional 
dependence of photon-initiated horizontal show- 
ers as opposed to neutrino-initiated showers that 
bear no evidence of the magnetic field |8I9| . 

The aim of this paper is primarily to discuss the 
experimental justification for the phenomenon 
and secondly to address the question of 'strong 
field' shower development in the many different 
magnetic fields which may be encountered by a 
photon of extragalactic origin. 



2. Pair production in a magnetic field 

In the following, pair production and radiation 
emission processes are considered as taking place 
in magnetic or electric fields that are homogenous 
(constant in space over the formation length of 
the corresponding process). The pair production 
probability differential in the energy of one of the 
final state particles, e, is given as ^01 
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where e/ = hu> — e, £ = 8u/3\, u = j 2 , and 
dW = dw/dt is the probability per unit time. 
Here e denotes the energy of one of the produced 
particles, e/ that of the other and Twj the energy 
of the photon. The functions if 1/3 and -K2/3 are 
the modified Bessel functions of order 1/3 and 
2/3, respectively. The Lorentz invariant strong 
field parameter is given as x — jB/Bq where 
B = m 2 c 3 /eh = 4.41 • 10 9 T is the critical 
magnetic field (corresponding to the electric field 
£0 = 1.32 • 10 16 V/cm) and 7 is understood here 
as hk±/mc, i.e. the photon momentum transverse 
to the magnetic field, in units of mc. The 'strong 
field regime' is reached when \ becomes compara- 
ble to or exceeds 1 - but already at x ~ 0.1 strong 
field effects become significant. Thus, for suffi- 
ciently energetic photons even the intergalactic 
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fields will appear as a strong field and generally 
speaking the Universe thus will become 'opaque' 
to these photons. 

The differential probability develops a pro- 
nounced minimum around x — e/huj = 1/2 and 
peaks at x ~ 1-6/x for large x- This - combined 
with the fact that the strong field radiation emis- 
sion tends toward the endpoint of the spectrum 
- means that as the energy increases beyond the 
region where \ is °f the order 100, the shower 
develops as essentially one energetic particle or 
photon, followed by many low energy particles. 
This behaviour is reminiscent of the behaviour 
of the Landau-Pomeranchuk-Migdal (LPM) effect 
for pair production (see e.g. |llj ) which also yields 
a preference for highly asymmetric pairs. 

The total pair production probability per unit 
time is given as JH] 

qm 2 c 4 f°° 8u + l 

W = TWT* / 3/2 / f g 2/3(fld« (2) 

6v37m w Ji u a ' A y/u—l 

In the limit x ^ 1 the probability is exponentially 
small, W ~ 3V3am 2 c 4 x/16V2n 2 wexp(8/3x), 
and in the limit % 3> 1 the result is W ~ 
O^amc^/Bo^X 1 ^ 3 , i- e - the probability actually 
diminishes once the energy surpasses the domain 
where x — 10 JUj. Indications of the 'satura- 
tion' precursory to the diminishment of the cor- 
responding energy loss by radiation emission has 
been shown experimentally only recently, |12I13| . 

The equations QJ and J2J are exact expressions 
(for B <C So), but eq. @ has an 'empirical ap- 
proximation' ^1] which generally forms the basis 
for the above mentioned calculations of the pair 
production in the Earth magnetic field. The de- 
cisive parameter for the behaviour is x = jB/Bo 
which depends only on the product of the field 
and the 'Lorentz factor' of the photon, Hk±/mc. 
The same x applies for strong field effects in ra- 
diation emission where x — E/mc 2 . 

With the acceleration techniques at hand to- 
day, the energy of electrons available for labo- 
ratory experiments is limited to 7 < 10 6 . This 
means that macroscopic magnetic fields of the 
order a few thousand Tesla or electric fields of 
the order 10 11 V/cm are needed to experimen- 
tally test the theory. At present this is certainly 
beyond reach by conventional means. However, 



an analogous case is seen for the pair production 
and radiation emission in single crystals and these 
phenomena enable a verification of the essential 
expressions, eqs. (QJ and 

3. Strong fields in crystals 

In the following, a brief introduction to strong 
field effects in crystals is given with particular em- 
phasis on pair production. However, essentially 
all conclusions and expressions can be applied 
also to the case of radiation emission and thus ap- 
ply for the development of electromagnetic show- 
ers as well. For a more complete introduction to 
strong field effects in crystals, see ^H], |T?S] 
and [T7j. 

3.1. Theory 

The large fields present near the nuclei in solids 
may in the case of single crystals add coherently 
such that a penetrating charged particle experi- 
ences a continuous field along its direction of mo- 
tion, see figure^ If further the particle is incident 
with a sufficiently small angle to a particular crys- 
tallographic direction, inside the so-called Lind- 
hard angle, the negatively/positively charged par- 
ticle is constrained to move near to/far from the 
nuclei and the electron clouds surrounding these. 
This is the channeling phenomenon |18| which has 
found widespread applications in physics. The 
critical angle for axial channeling is given as 

V>i = \/AZe 2 /pvd (3) 

where Ze is the charge of the lattice nuclei and 
d the spacing of atoms along the axial direction. 
The critical angle, V'l, is of the order a few tens 
of /xrad for e.g. 100 GeV electrons. Another well- 
known example of coherent action is the emission 
of so-called coherent bremsstrahlung. This ap- 
pears when crystal planes are crossed at regular 
intervals and at small angles to crystallographic 
directions during penetration [HJ\. The small an- 
gle requirement is essential for the coherence of 
the scattering that results in amplification of e.g. 
the radiation spectra. 

On the other hand, at high energies (generally 
above a few GeV for electrons) a seemingly dif- 
ferent coherent effect starts to appear. This ef- 
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of a few cm of a bent crystal. In this case, the 
equivalent field is as high as a few thousand Tesla. 

Strong field effects can be investigated by other 
means. One example is in heavy ion collisions 
where the field becomes comparable to the critical 
field, but the collision is of extremely short dura- 
tion. Another - technically demanding - example 
is in multi-GeV electron collisions with terawatt 
laser pulses where non-linear Compton scattering 
and so-called Breit- Wheeler pair production are 
observed [22] ■ In none of these cases can the field 
be considered macroscopic. 



Figure 1. A schematical drawing of the discrete 
nature of the scattering centers in a crystal and 
the resulting continuum approximation. 



feet arises due to the continuous field seen by 
the penetrating particle which acts on an elec- 
tron with sufficient strength to impose a strongly 
curved trajectory (as compared to the energy - 
and thus typical emission angles - of the parti- 
cle), see figure n The curvature - independent of 
the origin being magnetic or electric fields - gives 
rise to radiation emission of synchrotron charac- 
ter. Furthermore, this effect persists out to the 
so-called Baier angle, O , which is of the order of 
a few mrad and independent of energy 

6 ^ Ze 2 /dmc 2 (4) 

The electric fields, £, can be locally up to a 
few 10 11 V/cm, depending on the crystal type 
and orientation. The incident particle moves in 
these immensely strong fields over distances up to 
that of the crystal thickness, i.e. up to several cm. 
For a sufficiently energetic particle the critical pa- 
rameter x — t£/£q can thus reach values near 
and even beyond unity. Thereby the behaviour of 
charged particles in macroscopic strong fields as 
So or Bo can be investigated in accelerator based 
experiments. 

It is essentially the same phenomenon which 
is responsible for the mrad deflection of multi- 
GeV protons [20] and ions [2J during the passage 



3.2. Experiment 

In figurc[2]is shown the enhancement (the ratio 
of radiation lengths for the amorphous and crys- 
talline material) in a germanium crystal of the 
pair production probability differential in the en- 
ergy of one of the produced particles. This is com- 
pared with a calculation performed on the basis of 
eq. for mrad as shown by the line. The cal- 
culation is done by assuming a constant field dur- 
ing the formation of the pair and only an average 
over the fields encountered has to be performed. 
For the case of incidence with an angle of 2 mrad 
to the axis, the line denotes a strong field modi- 
fied coherent bremsstrahlung calculation |23| . In 
both cases a good agreement between experiment 
and theory is seen, details can be found in |24| . 

In figure 03 is shown the total enhancement in 
a tungsten crystal. The relevant averaged val- 
ues of x range up to about 5, i.e. strong field 
effects are expected to be dominant. The theo- 
retical values are based on eq. J2J) for the case on 
axis, whereas for higher angles the lines denote 
a strong field modified coherent bremsstrahlung 
calculation [23] • Again a good agreement be- 
tween experiment and theory is seen, details can 
be found in [2l)| . 

The slight disagreement at high photon ener- 
gies seen in figureEJis possibly due to the Landau- 
Pomeranchuk-Migdal (LPM) effect for pair cre- 
ation |27| which is not included in the calculation. 
The LPM effect is not noticeable until photon en- 
ergies of a few TeV in amorphous materials, but 
in crystals the multiple Coulomb scattering may 
increase by orders of magnitude and thus render 
the LPM effect visible at much lower energies. 
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Figure 2. Differential enhancement of pair pro- 
duction for photons in the energy interval 80-100 
GeV incident around the (110) axis in a germa- 
nium crystal. Filled squares are experimental re- 
sults and the line represents the theoretical val- 
ues. The results are shown for incidence in the 
(110) plane at angles and 2 mrad to the (110) 
axis 
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Figure 3. Totai enhancement of pair production 
for photons in the energy interval 10-150 GeV in- 
cident around the (100) axis in a tungsten crys- 
tal. Filled circles are experimental results and the 
dashed line represents the theoretical values. The 
results are shown for four angles to the (100) axis 
as indicated \2(tf . 
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Figure 4. The number of charged particles exit- 
ing an aligned 25 mm thick germanium crystal, 
compared to a Monte Carlo calculation (continu- 
ous line) based on eq. Q) and its equivalent for 
radiation emission \29lj . 



In figureQ]is shown the number of charged par- 
ticles exiting an aligned 25 mm thick germanium 
crystal. A comparison with a Monte Carlo calcu- 
lation based on eq. Q shows some disagree- 
ment, but in general the trend is correct. For 
details, see [2SI - It should be noted, however, 
that for the calculation of radiation emission in 
crystals a few additional complications arise due 
to a redistribution of the particle flux inside the 
crystal and the connected change in e.g. multiple 
Coulomb scattering. It is clear that small inaccu- 
racies in the assessment of the radiation probabil- 
ity or the pair production probability will prop- 
agate into amplified discrepancies after a large 
number of generations in the shower as could be 
the case for EASs. 

4. Pair production in the geomagnetic field 

In figure |S] is shown a calculation based on eq. 
similar to that of McBreen and Lambert [3| 
or Stanev and Vankov 0] to determine the con- 
version probability of a photon in the geomag- 
netic field. A dipolc field has been assumed with 
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Figure 5. The conversion probability of a photon 
in the geomagnetic Geld as a function of photon 
energy for a) energies in the region 10 1 — 10 3 EeV 
and b) energies in the region 10 5 — 10 17 EeV. 



values of 0.25 G and 0.528 G at the surface of 
the Earth, corresponding to the two sites of the 
Pierre Auger Observatory PHj. Figure repro- 
duces essentially the results of j.'illj for different 
fields while figure OJ) has been included to show 
the disappearance of the absorption for sufficently 
high values of \. The energy values used in fig- 
ure^ are certainly way beyond what is presently 
considered reasonable (even beyond the Planck 
energy) and serve only as an illustration of the 
W oc l/x 1 ^ 3 behaviour which may be relevant in 
other cases as discussed in the following. 



5. Pair production in extraterrestrial fields 

Following Hillas &\ there is a suggestive con- 
nection between the magnitude of extraterrestrial 
magnetic fields and their extent. Thus a simplifi- 
cation of the realistic fields can be parametrized 
roughly by 



B 



C 



(5) 



where C is set here somewhat conservatively to 1, 
but can be set to anything within 0.1 < C < 10 6 
- in any case the extent and magnitude of many 
extraterrestrial fields are uncertain [22]. The ex- 
tent of the field in parsec is given by L pc and B^q 
is the field strength in microgauss. By use of eq. 

the survival probability for a photon travers- 
ing regions with magnetic fields can be estimated. 
The result is shown in figure A change in C 
hardly affects the onset of the absorption band 
which is why such a large range of accepted val- 
ues can be replaced by a constant, set to 1 here. 
The explanation is that the decisive parameter 
is proportional to huB with an exponential in- 
crease, whereas the extent enters linearly. On the 
other hand, different values of C move the 'sur- 
vival island' in the upper, right corner since in 
this case there are two 'competing' power-laws. 
Thus, for relatively strong fields of limited extent, 
high energy photons may actually become more 
penetrating than those of lower energy. 

The range of figure has been set to include 
the possible value of an 'all-pervading' primordial 
field of < 10~ 9 G jS2] (where clearly the extent 
is severely underestimated) and to go near the 
highest field strengths at pulsars of ~ 10 s — 10 9 
T (but respecting B < B ). 

The above is certainly a rough simplification of 
the actual extraterrestrial fields which are in re- 
ality (as also seen from the 'Hillas plot') grouped 
in small clusters instead of being constrained to 
the line defined by eq. JSJ. This means that the 
'absorption band' shown in figure HJ is an over- 
estimate and should be replaced by 'islands' in- 
stead. However, this requires a thorough anal- 
ysis of the extent and magnitude of the actual 
fields and is beyond the scope of this paper. Fur- 
thermore, there are other absorption mechanisms 
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Figure 6. The conversion probability as a func- 
tion of photon energy and magnitude of the mag- 
netic field under the assumption implicit in eq. 
© with C = l. 



for photons: Collisions with the infrared/optical, 
cosmic microwave background or universal radio 
background photons [5]. 

6. Conclusion 

It has been shown in the presented paper that 
the theoretical approaches used for simulations of 
photon cascade development in e.g. the geomag- 
netic field are experimentally justified from the 
point of view of 'strong field effects' in crystals. 

Actual experiments using crystals may to a 
large extent simulate the behaviour of shower 
development in the Earth magnetic field and 
atmosphere as a number of effects (strong 
field, LPM, Chudakov as well as e.g. ordinary 
bremsstrahlung) are believed to be present in 
both cases. In fact, model builders may use data 
obtained in crystals as test cases for the simula- 
tion of many aspects of extended air showers. 
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